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Abstract 

In this paper, we investigate the use of relaxed logarithmic barrier functions in the context of linear model predictive 
control. We present results that allow to guarantee asymptotic stability of the corresponding closed-loop system, and 
discuss further properties like performance and constraint satisfaction in dependence of the underlying relaxation. The 
proposed stabilizing MFC schemes are not necessarily based on an explicit terminal set or state constraint and allow 
to characterize the stabilizing control input sequence as the minimizer of a globally defined, continuously differentiable, 
and strongly convex function. The results are illustrated by means of a numerical example. 


I. Introduction 

T He basic idea of model predictive control (MPC) is to 
obtain a stabilizing feedback law by solving at each 
sampling instant, in a receding horizon fashion, a suitable 
finite-horizon open-loop optimal control problem that is 
parametrized by the evolving system state. Only the first 
element of the optimal control input is applied to the 
plant and the optimization is repeated at the next time 
step in a receding horizon fashion. Based on the pre¬ 
dicted dynamics of the system to be controlled, both a 
user-defined cost objective and potential constraints on 
the system states and input can be taken into account 
explicitly within the corresponding optimization problem. 
There exist various theoretical results concerning the sta¬ 
bility properties of the closed-loop system for both linear 
and nonlinear systems, and MPC can be considered as a 
widely accepted control concept that is more and more 
applied to industrial processes, see 221 and [i^ as well as 
references therein. 

However, there often is a gap between theoretical MPC 
concepts and their actual implementation, and several 
disadvantages remain. For example, algorithmic aspects 
of the required optimization are usually not taken into 
account explicitly in the MPC design, which may result 
in suboptimal or even unstable closed-loop behavior. Fur¬ 
thermore, in the presence of disturbances, sensor outliers, 
or observer errors, the underlying open-loop optimal con- 
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trol problem may become infeasible, leading to a com¬ 
plete crash of the respective control algorithm. In recent 
years, considerable research effort has been spent on de¬ 
signing efficient and reliable algorithmic MPC implemen¬ 
tations that, if possible, still preserve some of the desired 
theoretical properties when running in closed-loop oper¬ 
ation, see e.g. iii [3^ . One particular approach 
that allows to reduce the gap between MPC schemes on 
the one and MPC algorithms on the other hand is given 
by the concept of barrier function based MPC, which 
was introduced in and has been extended recently 
in [l3, [m, [ 3 . In barrier function based MPC schemes, 
the inequality constraints occurring in the open-loop op¬ 
timal control problem are incorporated into the cost func¬ 
tion by means of suitable barrier function terms like it is 
also done in interior-point optimization algorithms 241. 
As shown in the above references, barrier function based 
MPC approaches allow to guarantee asymptotic stability 
of the closed-loop system while at the same time reduc¬ 
ing the underlying open-loop optimal control problem to 
an equality constrained or even unconstrained optimiza¬ 
tion problem. Thus, a barrier function based reformu¬ 
lation is already integrated into the MPC design, which 
then allows for an efficient implementation based on tai¬ 
lored numerical optimization techniques. In particular, 
the barrier function based approach makes the open-loop 
optimal control problem accessible for continuous-time al¬ 
gorithms that asymptotically track the optimal solution 
of the unconstrained reformulation and, thus, allow to 
implement MPC without any iterative on-line optimiza¬ 
tion, see However, two main disadvantages 
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remain. On the one hand, the open-loop optimal con¬ 
trol problem may still become infeasible since the barrier 
functions are only defined within the interior of the corre¬ 
sponding constraint sets. On the other hand, the barrier 
function based approach inherently requires the use of a 
terminal set or terminal equality constraint, which is not 
desirable, and also typically not used, in practice. 


In this paper, we show that both of these problems can 
be eliminated by making use of so-called relaxed loga¬ 
rithmic barrier functions, i.e., barrier functions that are 
smoothly extended by a suitable Mnalty term outside of 
the corresponding constraint set [1, [13,113 ■ We discuss 
suitable relaxation procedures (Section]^ and present re¬ 
sults that allow to guarantee asymptotic stability of the 
closed-loop system with and without making use of ter¬ 
minal sets fSectionsiBl ICj) . to determine and control the 
maximal violation of input and state constraints in closed- 
loop operation (Section [D|, and to recover the optimal¬ 
ity properties of both nonrelaxed barrier function based 
and conventional linear MPC schemes (Section [Q. Fur¬ 
thermore, based on the presented results, we provide in 
Section El a step-by-step procedure for the constructive 
design of the overall MPC scheme. 

The basic idea of relaxed barrier function based MPC is 
closely related to approaches based on soft constraints 
or penalty functions, see e.g. [13, js^l- However, in these 
approaches, closed-loop properties like stability or strict 
constraint satisfaction can usually only be guaranteed 
when making use of nonsmooth or even exact penalty 
functions. The key feature of the relaxed barrier functions 
discussed in this paper is that they result in a smooth for¬ 
mulation of the overall problem while still allowing for an 
arbitrary close approximation of the corresponding non¬ 
relaxed case. As a result, the stabilizing control input 
can be characterized as the minimizer of a globally de¬ 
fined, twice continuously differentiable, and strongly con¬ 
vex cost function, which makes it efhciently computable 
by standard nonlinear programming algorithms. Some 
preliminary results on the concept of relaxed barrier func¬ 
tion based MPC have been presented in [l3|. In this pa¬ 
per we significantly extend these results and provide an 
in-depth study of several theoretical and practical aspects 
of the presented linear MPC schemes. A particularly in¬ 
teresting new result is the insight that the use of relaxed 
barrier functions allows to ensure global asymptotic sta¬ 
bility of the closed-loop system without making use of a 
terminal set. 


Throughout the paper we will make use of the follow¬ 
ing notation. R_|_, K+|-, and N+ denote the sets of non¬ 
negative real, strictly positive real, and strictly positive 
natural numbers. Furthermore, §" and refer to the 
sets of positive semi-definite and positive definite matri¬ 
ces of dimension n G N+. For any given matrix M or 
vector V, and u® refer to the i-th row or element and 
\\x\\m := Vx^Mx for any M £ S+. A polytope is defined 
as the compact intersection of a finite number of halfs¬ 


paces. For any arbitrary set S, the expression S° will de¬ 
note the open interior and dS the boundary. Moreover, 
I := [l • • • l] of suitable dimension. 

II. Problem Setup 

In this paper, we consider the control of linear time- 
invariant discrete-time systems of the form 

x{k + I) = Ax{k) + Bu{k ), (1) 

where x{k) £ R” refers to the vector of system states and 
u{k) £ R”^ refers to the vector of system inputs, both at 
time instant k > 0. Moreover, the matrices A £ R”^" 
and B £ describe the corresponding system dy¬ 

namics, where we assume (A, B) to be stabilizable. The 
control task is to regulate the system state to the origin 
while minimizing a given, user defined performance crite¬ 
rion and satisfying state and input constraints of the form 
x{k) £ X and u{k) £ U for all A: > 0. Here, X C R” and 
U C R"* are typically considered to be given convex sets 
that contain the origin in their interior. In this paper, we 
assume X and U to be polytopes defined as 

A = {a; £ R" : C^x < dx}, (2a) 

U — {uG R"* : CuM < du} , (2b) 

where Cx £ R?-''", £ R9“X'" and dx £ R^, du £ R^ 

with qx,qu G N+. In linear MPC, this problem setup is 
usually handled by solving at each sampling instant an 


open-loop optimal control problem of the form 

N-l 

/ ^{^k, Uk) + F{xn) (3a) 

s.t. Xk+1 = Axk + Buk, xq = X , (3b) 

Xk G X, k = 0, ..., N — 1, xn G Xf , (3c) 

Uk Gl4, k = 0,.. ., N — 1 . (3d) 


for the current system state x = x{k) and a finite predic¬ 
tion and control horizon N G N+. Here, the stage cost 
£ : R” X R™ —^ R-|_ and the terminal cost F : R” —>• R+ 
are usually defined as i{x,u) = + ||u||^ and F{x) = 

||a;||p for appropriately chosen weight matrices Q £ S”, 
R G P G S!)^. Furthermore, u = {uq, ..., mat-i} 
denotes the sequence of control inputs over the predic¬ 
tion horizon N, while Xf refers to a closed and con¬ 
vex terminal constraint set that may be used to guar¬ 
antee stability properties of the closed-loop system. Note 
that we make use of subindices to distinguish open-loop 
predictions Xk, Uk from actual state and input trajecto¬ 
ries x{k), u{k). The control law is obtained by solving ([3]) 
at each sampling instant A: > 0 and applying u{k) = 
UQ{x{k)) in a receding horizon fashion. Sufhcient con¬ 
ditions for the recursive feasibility of ([3]) as well as for the 
asymptotic stability of the closed-loop system are summa¬ 
rized in [ 2 ^. In the following, we will refer to this setup 
as conventional linear MPC. Moreover, we will often use 
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to denote the successor state and write © as x'^ = 
Ax + Bu. For given system state x = x(k) and input se¬ 
quence u = {uq, ..., itAT-i}, the resulting open-loop state 
sequence is given by x{u,x) = {xq{u,x), ... ,xn{u,x)'\, 
where the elements Xk{u,x), k = 0,... ,N are given ac¬ 
cording to (ISb) with xq{u,x) = X. For a given optimal 
input sequence u*(x), we write x*{x) := x{u*(x),x) as 
well as xl(x) := Xk(u*(x), x). Sometimes we also drop 
the explicit dependence on the current system state for 
ease of notation. 


III. Preliminary Results 

We introduce in this section some concepts and results on 
nonrelaxed logarithmic barrier function based MPC that 
we will need in the remainder of the paper. In particular, 
we present a novel stability theorem which generalizes 
different existing ideas and is applicable to a wide class 
of barrier function based MPC approaches. 


A. Barrier Function Based Model Predictive Control 


The main idea in barrier function based MPC is to elimi¬ 
nate the inequality constraints from the above MPC open- 
loop optimal control problem by making use of suitable 
barrier functions with a corresponding weighting factor. 
Based on this idea, it is possible to reformulate prob¬ 
lem ([3]) as an equality constrained (or even unconstrained) 
strongly convex optimization problem, which can then be 
solved by means of tailored optimization procedures like, 
e.g., the Newton-method. In general, the exact solution 
to the original problem, and hence also the stability prop¬ 
erties of the corresponding closed-loop system, are recov¬ 
ered when the weighting factor of the barrier functions 
approaches zero [a, ch. 11]. However, for an arbitrary 
but fixed nonzero weighting, as it will necessarily occur 
in all numerical implementations, stability of the origin is 
by no means guaranteed. Several approaches towards the 
stabilizing design of barrier function based MPC schemes 
have been presented in [3, 11, [^. In the follow¬ 


ing, we shortly summarize the main aspects of barrier 
function based linear MPC and present a fairly general 
stability theorem for the considered problem setup with 
polytopic input and state constraints. Details on the dis¬ 
cussed MPC approaches can be found in the respective 
references. 

Let us consider in the following the barrier function based 
open-loop optimal control problem 


N-l 

Jn{^) = min i{xk^ Uk) + F{xn) (4a) 

s. t. Xk+1 = Axk + Buk, Xo = X , (4b) 

where £{x,u) := £{x,u) + eB^{u) + eB^{x) and F{x) := 
F{x) + eB{{x) are the modified stage and terminal cost 
terms for £{x,u) = HxHq -I- ||u|||j and F{x) = ||a;||p, and 


-Bu(‘)) ^x(‘)) ^nd B[{-) are suitable convex barrier func¬ 
tions with domains U°, X°, and Xj with Bu{u) —>■ oo 
for u —>■ dU, B^{x) —>■ oo for X —>■ dX, and Bf{x) —>■ oo 
for X —^ dXf. Furthermore, the positive scalar e > 0 
is the barrier function weighting parameter which deter¬ 
mines the influence of the barrier function values on the 
overall cost function. As outlined above, the goal is now 
to choose the problem parameters P and Xf a.s well as 
convex barrier functions i?u('): ^x(’): ^^.d B[{-) in such 
a way that a linear MPC scheme based on asymp¬ 
totically stabilizes the origin for any arbitrary but fixed 
e S K.-I+. If we want to employ standard MPC stabil¬ 
ity concepts which are based on using the value function 
J)(r(-) as a Lyapunov function for the closed-loop system, 
this usually requires the overall cost function to be con¬ 
tinuous as well as positive definite with respect to the 
origin. In order to ensure the latter property for general 
barrier functions, the concept of gradient recentered bar¬ 
rier functions has been introduced in Sg- In addition, a 


weighting based recentering approach for logarithmic bar¬ 
rier functions on polytopic sets has been proposed in |13|. 


Definition 1 (Recentered log barrier function [iSlIl). 
Let V = {z € ML : Cz < d} with C S d € R.^ 

be a polytopic set containing the origin and let B :V° ^ 
R, B{z) = Bi{z) with Bi{z) = — ln(—-I-d®) he 
the corresponding logarithmic barrier function. Then, the 
function Bq ■ V° —>■ R+ with 


Bg{z)=B{z)-B(Q)-[VB{Q)]^z 


( 5 ) 


defines the gradient recentered logarithmic barrier func¬ 
tion for the polytopic setP. Furthermore, Bw ■ P° R-i- 
with 

9 

Rh^(z) =^(l + u;*)(R,(z)-R,(0)) (6) 

i=l 

defines the corresponding weight recentered logarithmic 
barrier function, where the weighting vector w € R^ 
is chosen in such a way that VBw{0) = + 

w^)VB,{0) = 0 . 

In principle, the recentering can be seen as a modi¬ 
fication which preserves the main characteristics of the 
barrier function while ensuring that it is positive definite 
with respect to the origin. Both recentering approaches 
can in general be applied to our problem setup and will 
result in a continuously differentiable cost function that is 
positive definite and strongly convex in the optimization 
variable u. Moreover, the stability results presented in 
this work do hold independently of the underlying recen¬ 
tering method. 

Let us now turn towards the stability properties of the 
closed-loop system when applying the barrier function 
based MPC feedback u{k) = u{x{k)). 

Definition 2. For N G N+, let the feasible set Xm he 
defined as Xn := {x G X : 3u = {uq, ...,uat-i} such 
that Uk G U, Xk{u, x) G X for k = 0,..., IV — I and 
xn{u,x) G Xf}. 
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Definition 3. In the following, the matrix A k ■= A+BK 
describes the closed-loop dynamics for a given stabilizing 
local control law u = Kx and Bk{x) := Bx{x) + Bn{Kx) 
refers to the corresponding combined barrier function of 
input and state constraints for the set Xk := {x £ X : 
Kx G U}. 


As outlined above, different approaches towards the 
stabilizing design of barrier function based MPC formu¬ 
lations, i.e., on how to choose the terminal cost matrix P, 
the terminal set Xf, and the correspon ding barrier func¬ 
tion B{{-), have been presented in and [l^ . 

In principle, all mentioned approaches are based on the 
idea of choosing the terminal set as a positively in¬ 
variant subset of the state space in which the function 
Bx{-), i-C., the influence of the input and state constraint 
barrier functions, can be upper bounded by a quadratic 
function. Then, based on this quadratic bound, the ter¬ 
minal cost matrix P is computed in such a way that it 
compensates for this influence and ensures that the bar¬ 
rier function based terminal cost F{-) is a local control 
Lyapunov function for the auxiliar y c ontrol law u = Kx. 
While the approaches in and [1^ make use of ellip¬ 
soidal terminal sets, the approaches presented in [HI and 
[l^ are based on poly topic terminal sets. In the follow¬ 
ing, we present a set of sufficient stability conditions that 
is used in the remainder of this work and can be seen 
as a generalization of the ideas presented in the above 
references. 


Assumption 1. For Q G i? G e G R.-|+ and 

a given stabilizing local control gain K G let the 

parameters of the barrier function based open-loop optimal 
control problem satisfy the following conditions. 


Al: The barrier functions B^(-) and i?x(‘) cire recentered 
barrier functions for the sets U and X, respectively. 

A2: There exists M G s. t. Bk{x) < x^Mx Mx G N, 
where M C Xk, 0 G M°, is a convex, compact set. 

A3: The matrix P G is a solution to the Lyapunov 
equation P = Aj^PAj^ -\- K^RK -|- Q -b sM. 

A4: The terminal set Xf is convex and compact with 0 G 
Xf, Xf C Af C Xk, and = Akx G A/ Vx G Xf. 

A5: The function B{{-) is a recentered barrier function 
for the set Xf and B{{Akx) — Bf{x) <0 Wx G Xj. 

Based on Assumption [T] we can state the following re¬ 
sult on the stability of the closed-loop system. 


Theorem 1. Let Assumption [7] hold true and let the 
feasible set Xk defined according to Definition [H have 
nonempty interior. Then, the barrier function based 
MPC feedback u{k) = UQ{x{k)) resulting from problem m 
asymptotically stabilizes the origin of system m under 
strict satisfaction of all input and state constraints for 
any initial condition x(0) G Xfj. 


Proof. The proof uses standard MPC stability arguments 
and comprises some of the main ideas presented in [3113 


and mi. We first show recursive feasibility of prob¬ 
lem 0. For any Xq G Xf^, there exists by definition 


an optimal input sequence u*(xo) = {u. 


that 


0) • ■ • ) “-AT 

guarantees strict satisfaction of all input, state, and ter¬ 
minal set constraints and results in a feasible open-loop 
state sequence x*(xo) = {xo,x*,... ,x^} with Xk G X^. 
The successor state is given as Xq = Axq -\- Buq = x^. 
Due to the linear dynamics and the properties of the 
terminal set Xf, see A4 of Assumption [U ■u^(xo) = 
{mj, ... ,u*K_i, Kx*k} is a suboptimal but feasible input 
sequence for the initial state Xq that results in the in 
the feasible open-loop state sequence a;(it'''(xo),x[}") = 
{xq ,X 2 ,... ,x*j,f, AKx"ff} with Akx'^ G XJ. This shows 
that for any xq G Xff, there exists a feasible input se¬ 
quence that ensures that the successor state Xq = Axq -b 
Buq lies again in the interior of the feasible set Xk, which 
guarantees recursive feasibility of the open-loop optimal 
control problem (|4]). 

In the following, we show that the value function satisfies 


J*Ni4) - J*Nixo) < -i{xo,Kixo)) Vxo G . (7) 

First, due to the suboptimality of u'^{xo), it holds that 
Jjfix^) - Jk{xo) < JAr(w’'‘(xo),xJ) - Jk^xq), where 
Jk{u~^{xo),Xq) denotes the value of the cost function 
evaluated for the suboptimal input sequence 'U~'’(xo). 
Moreover,^ Jjv(m+(xo),x,|) - Jf^ixo) = F{Akx%) - 
F{x*j,^) -b i{x%,Kx*j.f) - l(xo,'iio) < -i{xo,Uo) for any 
xo G Xfj since 

F{Akx*k) - F(x^) + i{x*K, Kx*k) (8a) 

= \\Akx*k\\1 - \\x%\\l + Ik^vllg + (8b) 

-b ei?if (xjv)-b ei?f (Aa'Xjv) — ei3f(x)v) (8c) 

< e {Bi{Akx*k) - Bi{x*K)) <0 ^ x*k G X] . (8d) 


Here, the first inequality follows from the quadratic bound 
Bk{x) < 'yx^Mx yx G Xf c Af and the suitable choice 
of the terminal cost matrix P, see A2 and A3 of Assump- 
tion[T] The second inequality holds due to A5 of Assump¬ 
tion [TJ 

Finally, the considered problem setup and the design of 
the barrier functions according to Al of Assumption [T] 
ensure that : X^ —>• R+ is a well-defined and positive 
definite function with J)(r(x) —>■ oo whenever x —>■ OXk. 
Thus, in combination with the decrease property 0, it 
can be used as a Lyapunov function which proves asymp¬ 
totic stability of the origin of system o under the feed¬ 
back u{k) = UQ{x{k)) for any x(0) G Xff. □ 

Remark 1. Note that the barrier function parameter 
e G K 4 f can in principle be decreased iteratively in each 
sampling step or between two consecutive sampling steps 
without loosing feasibility or stability properties of the 
closed-loop system. This might be meaningful for numer¬ 
ical reasons or in order to enforce convergence to the op¬ 
timal solution of the conventional MPC problem (0)- We 
limit ourselves to fixed values of e. 
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Different approaches on how to actually construct the 
neighborhood M and the corresponding quadratic bound 
for Bk{‘) as well as suitable choices for the terminal set 
Xf and the barrier function -Bf(-) have been proposed 
in [SElIiSlal. However, the fact that the underlying 
barrier functions are only defined in the interior of the re¬ 
spective constraint sets may be problematic both from a 
practical and from a conceptual point of view. On the one 
hand, violations of the state, input, and terminal set con¬ 
straints are not tolerated at all, which might cause severe 
problems in the presence of uncertainties, disturbances, 
noise, observer errors, or sensor outliers. On the other 
hand, all existing stability concepts inherently require the 
use of a suitable terminal set as they are based on upper 
bounding the barrier function Bk{-) by a quadratic func¬ 
tion, which is of course only possible locally in a region 
around the origin. In the following section, we introduce 
the concept of relaxed logarithmic barrier function based 
MFC and show that it can be used to overcome all these 
limitations, allowing for conceptually simpler and more 
reliable linear MFC schemes. 


IV. Main Results 


The basic idea of relaxed logarithmic barrier functions is 
to smoothly extend a given logarithmic barrier function 
with a suitable, globally defined penalty function 
[ 2 ^ . In the following, we provide an in-depth study of sev¬ 
eral interesting theoretical and practical aspects of linear 
MFC approaches that are based on such relaxed logarith¬ 
mic barrier functions. In particular, we show that, on 
the one hand, feasibility and stability properties of the 
nonrelaxed formulation can always be recovered by ap¬ 
proximating the original barrier functions close enough. 
On the other hand, we present novel linear MFC schemes 
that, by exploiting the properties of relaxed logarithmic 
barrier functions, do not require the use of a terminal set 
and allow to guarantee global asymptotic stability of the 
closed-loop system, see Section [Cl The benefits of global 
stabilization and a simplified design procedure are bought 
with the loss of guaranteed input and state constraint 
satisfaction. However, as we will show in Section El the 
relaxed barrier functions can always be designed in such 
a way that an a-priori defined tolerance for the maximal 
violation of state and input constraints is guaranteed for 
a certain set of initial conditions. We briefly discuss some 
interesting performance and robustness properties of the 
proposed MFC schemes in Section El 


A. Relaxed Logarithmic Barrier Function Based MFC 

We begin our studies by introducing the concept of re¬ 
laxed logarithmic barrier functions and discussing suit¬ 
able realizations based on different relaxing functions. 

Definition 4 (Relaxed logarithmic barrier function). 
Consider a scalar 6 S R 4 f 7 called the relaxation param¬ 
eter, and let P{-;S) : (— 00 , <5] —>• R 6e a strictly mono¬ 


tone and continuously differentiable function that satisfies 
P{5] 5) = — ln(d) as well as lim 2 _>._oo l3{z; 6) = 00 . Then, 
we call R : M —>■ R defined as 


B{z) 


— ln(z) z > 6 

l3{z; 6) z <5 


(9) 


the relaxed logarithmic harrier function for the set R_|_ and 
refer to the function as the relaxing function. 

A graphical illustration of the basic idea is given in 
Fig. [TJ In general, it is advisable to choose the relax¬ 
ing function /?(• ;(5) as a strictly convex function that 
smoothly extends the natural logarithm at z = d. In 
this case, Bf) is a strictly convex function that is twice 
continuously differentiable and defined on z G (— 00 , 00 ). 
Note that lim 5 ^o-B(z) — ^ B{z) for any strictly feasible 
z G IR-h-, which shows that the nonrelaxed formulation 
can always be recovered by decreasing the relaxation pa¬ 
rameter 8 to zero. Note that we do not indicate the ex¬ 
plicit dependence of the relaxed barrier functions on the 
relaxation parameter 5 for the sake of notational simplic¬ 
ity. However, we will use Bf) to denote the relaxed ver¬ 
sion of a barrier function based expression R(-). 

The first ideas on relaxed (or approximate) logarith¬ 
mic barrier functions with a quadratic relaxing function 
/3(- ; 8 ) seem to have been proposed in i and [ 2 ^, re¬ 
spectively. In 0, the authors extended the concept to 
general polynomial penalty terms and applied it in the 
context of continuous-time trajectory optimization. In 
particular, the authors make use of the polynomial relax¬ 
ing function 


Pk{z-,8) 


k — 1 f z — k8 


ln(d), (10) 


where fc > 1 is an even integer. It is easy to verify that the 
function /3fc('; i5) has all the desired properties mentioned 
above. As reported in [17|, already a quadratic extension 
based on fc = 2 seems to work well in practice. 

In order to avoid large constraint violations, it may be 
beneficial if the relaxing function increases very rapidly 
outside the border of the feasible set. As an alternative 
to the polynomial relaxation above, we proposed in [l^ 
the following exponential relaxing function 


fie{z] 8) = exp (^1 - - 1 - ln((5), 


( 11 ) 


which is an upper bound for the function /3k(-',8). In 
fact, using the derivatives /3(,(-;d), /?e(-;8) and the limit 
representation of the exponential function, it can be 
shown that fik(z;8) < fik-i-2(z;8) < fie(z;8) and that 
limk->.oo fik(z; 8) = fie(z;8) Vz < 8. Furthermore, for 
/3(-; 8) = Pk {-; 8) and /3(- ;8) = fief 8) it also holds that 
B(z) < B(z) Vz G R^. We assume in the following that 
one of the two functions fikf ;d) and fief ;8) is used as 
relaxing function. 
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Based on these concepts, let us now consider the following 
relaxed barrier function based MFC formulation 


N-l 

J'^{x;S) = uimy^ i{xk,Uk) + F{xn) (12a) 

s.t. Xk +1 = Axk + Buk, Xo = X , (12b) 

where £{x,u) := £{x,u) + eBu{u) + eB^{x) for suitable re¬ 
laxed recentered logarithmic barrier functions i?u(’) and 
i?x(-) as discussed below. The term F{-) denotes a suit¬ 
able relaxed terminal cost function which we do not spec¬ 
ify at the moment. As we will see in the Sections|B]and O 
the choice of F{-) is crucial for ensuring stability proper¬ 
ties of the closed-loop system. 

Assumption 2. The functions B^j^{u) = J2i=i 
and B^{x) = 'TAAi^^AA relaxed recentered loga¬ 
rithmic barrier functions for the polytopic sets U and X 
and the relaxation parameter 5 G is chosen such that 
B^{0) = B^(0) = 0 and VB^iO) = 0, V.Bx(0) = 0. 

Note that a S satisfying Assumption [5] always exists. 
In particular, if the relaxation parameter satisfies 0 < 
6 < min{dO,..., d),, ..., d®“}, suitable relaxed barrier 

functions can be easily constructed by simply relaxing the 
logarithmic terms of the recentered nonrelaxed formula¬ 
tion. Hence, for Zi{x) := —C^x -|- d^, i = 1,..., Qx, we get 
for example 

- /N i-Hz^ix))+HdA - ^ z^(x)>s 

B.Ax)={ . (13) 

[/3{zi{x)]6) -\-\n{dD - ^ Zi{x)<5 

when using gradient recentered barrier functions and 

B f(l + 0 +M4)) Zi{x)>5 

\(1-bwi) (^(zi(a;);(5)-|-ln(d(,)) Zi{x) < 5 

(14) 

when using weight recentered barrier functions, cf. Defi¬ 
nition [TJ The barrier functions B^^i{-) for the input con¬ 
straints can be defined analogously. Note that it is in prin¬ 
ciple possible to consider individual relaxation parameters 
6i for the different constraints, but we restrict ourselves 
to one overall <5 for the sake of simplicity. While the relax¬ 
ation of the barrier functions directly implies that prob¬ 
lem (fT^ is always recursively feasible, stability properties 
of the resulting closed-loop system are by no means guar¬ 
anteed. This problem will be addressed in the following 
two sections. 


B. Closed-Loop Stability: Terminal Set Based Ap¬ 
proaches 


In this section we present our main results on the closed- 
loop stability properties of relaxed logarithmic barrier 
function based MFC schemes that make use of a suitable 
terminal set, cf. 12 1. In accordance with [i3,[Iil,[i3,[3i|, 
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Figure 1: Left: Principle of relaxed logarithmic barrier 
functions based on the relaxing function P{z;5). Right: 
Weight recentered logarithmic barrier function (solid) and re¬ 
laxed weight recentered logarithmic barrier function for S £ 
{0.01, 0.1, 0.5, 1} for the constraint —1 < 2 < 2 with 2 £ R. 


we assume in the following that the terminal set Xf can 
be represented as 


Xf = {x G : ip{x) < 1} , 


(15) 


where (p : K" —>■ K+ is a continuously differen¬ 
tiable, convex, and positive definite function that satis¬ 
fies p{Akx) < <p{x) Vx £ Xf. When considering ellip¬ 
soidal terminal sets as discussed in and[l^. pf) may 
for example be chosen as pix) = x^ PfX with a suitable 
Pf G In ll| and 13[, it is moreover shown how 

smooth approximations of the Minkowski functional can 
be used to dehne p(-) for polytopic terminal sets of the 
form Xf = {xGU'^: HfX < 1}. 


B1) Stabilization with strict constraint satisfaction 

In the following, we will show that asymptotic stability 
of the closed-loop system as well as strict satisfaction of 
all input and state constraints can always be guaranteed 
for any feasible initial condition by making the relaxation 
parameter S G arbitrarily small. The key assumption 
underlying the following results is that both the state and 
input constraint barrier functions and the terminal set 
constraint barrier function are relaxed. In particular, the 
barrier functions for the input and state constraints are 
chosen according to Assumptionj^l while the terminal cost 
function is given by 

F{x) = x^Px-\-eBf{x) , (16) 

where B{(-) is a relaxed logarithmic barrier function for 
the terminal set Xf. Assuming 5 < I, the function Bi{-) 
can be defined as follows. 

Assumption 3. The barrier function Bi{-) for the ter¬ 
minal set Xf is a relaxed logarithmic barrier function of 
the form 


B{{x) 


— ln(l — p{x)) 1 — p{x) > 5 

I3{l-p{x)]5) l-p{x)<5, 


(17) 


where 5 G 5 < 1, is the relaxation parameter and 

p{-) is the function defining the terminal set according 

to 1751) . 
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Note that the function Bf(-) is continuously differen¬ 
tiable, convex, and positive definite by design. Consider 
now the following definition, which introduces a lower 
bound for the value of the relaxed barrier functions 

and Bi{-) evaluated at the borders of the respective 
constraint sets. 

Definition 5. Let the scalar /j(S) € R+|- be defined as 

/3((5) = min{/3x(i5),/3u(^),/3f((5)} , (18) 

where = /3(0;d) and 

= inm{B^{x)\ C^x = d^}, i = (19a) 

i,x 

fiuiS) = mm{Bu{u)\ = d^}, j = l,...,qu- (19b) 

U,j 

Note that the values Pu{b), and hence also i3{S), 

can be computed easily for a given 6 G as the op¬ 
timization problems in m are convex. Moreover, for 
the case of gradient recentered relaxed logarithmic bar¬ 
rier functions, an explicit expression for a lower bound on 
fi{6) has been given in [l^ . The following Lemma shows 
that the sublevel sets of iix(-)! ^u(-)i Bf{-) related to 

fixib), /3u(<5)) and are always contained within the 

sets X, U, and Xf, respectively. 

Lemma 1. Let the Assumptions and\^ hold and let 
the values l3u{S), and /3f(5) be defined according to 

Definition\^ Then it holds that = {x G R”|.Bx(3^) < 

^x(5)} C T”, = {u€ R"^\bS) < ^u(5)} C u, as 

well as = {x G R"|.Bf(x) < ;5f((5)} C Xf. 

Proof. The result for the terminal set Xf follows directly 
from the fact that Bf{-) is convex and strictly monotone 
in the argument z = 1 — (p{x). Considering the state 
constraints, let us assume that there exists an x G Sb^ 
with X ^ T”, i.e., CxX > dx- However, then there would 
exists a A G (0,1) such that CxAx = XC^x < dx and 
C^Ax = XC^x = dx Ini' some i = 1, ..., gx- Now, due to 
the convexity and positive definiteness of Hx(-), it holds 
that Bx{Xx) < XBx{x) < which obviously is a 

contradiction to the definition of /3x(d). For the input 
constraints we can use similar arguments. □ 

Based on the above results, we now define a set of initial 
conditions for which we can guarantee asymptotic stabil¬ 
ity of the origin as well as strict satisfaction of all input 
and state constraints. 

Definition 6 . For a given 6 G and a corresponding 
j3(5) G K.-I+ according to Definitions^ let the set Xn{S) 
be defined as Xn{ 6) := {x G M" | Jff{x] 6) < efi{S)}. 

Theorem 2. Let Assumptions [H [H and [51 hold true. 
Consider problem m with the terminal cost F{-) 
from (OH) and let the set Xj\j{S) be defined according to 
Definitions^ Then, the feedback u{k) = UQ{x(k)) asymp¬ 
totically stabilizes the origin of system m under strict 
satisfaction of all input and state constraints for any 
x(0) G Xn{S). 


Proof. The proof consists of three parts and is closely 
related to that of Theorem 12 in 12|. First, we show 
that the underlying input, state, and terminal set con¬ 
straints are not violated for any xq G Xn{S); then we use 
standard MFC arguments to show that the value function 
J*{x{k);S) will decrease when applying the feedback u{k)-, 
finally, we use this result to conclude that the resulting 
input and state sequences will also be strictly feasible at 
all later time steps and that the origin of the closed-loop 
system is asymptotically stable. 

-il, x*{xo) = {xo,.. • ,x;^} 


-0, . . . , U.JV 


i) Let u*(xq) = {? 
denote the optimal open-loop input and state sequences 
for a given xq G XN{b). Since the cost function in (fT^ is 
a sum of positive definite terms, it holds that eBx{xl) < 

sBuSul) < Jff{xo]6), as well as ei?f(x^) < 
Jff{xo; S) for all xq G T’Ar((5) \ {0}. Hence, Bx{xl) < i3{S), 
Buiul) < fi{S), as well as Bfix'^) < fi{6). Due to 
Lemma |T] and the definition of /3(d), this implies that 
xl G X° and ul G U° as well as x*^ G Xj for any 
Xo G <Tjv(d) \ {0}. Of course, the case xq = 0 is triv¬ 
ial. Hence, the predicted input and state sequences are 
strictly feasible and the applied input results in a succes¬ 
sor state x)j" = Hxo -I- Bii^ G Xf^. 

ii) Let us now consider x[j" G Xfj. We can use basically 
the same arguments as in the proof of Theorem |T] to show 
that 


J)(f(x)j';d) - i^(xo;d)<-i(xo,Mo(a^o)) Vxq G TW(d). 

( 20 ) 

In particular, we know that it^(xo) = 
{u\,... and x+(xo) = 

{x/}", X 2 ,..., x)/^, with Akx*j.^ G Xf are sub- 

optimal but feasible input and state sequences for the 
initial state x/J" G X^. Moreover, 

FiAKX*j,f)-Fix%)+iixlf,Kx%) ( 21 a) 

= ||H;,x^fp - llx^lll, + ||x)Vll^ + \\Kx*j,rR ( 21 b) 
sBx{x%) sB{(Akx'^) — sBfSx*]^) ( 21 c) 

< e{MAKX%) - Mx*n)) <0 'i x*x€Xf . ( 21 d) 

Here, the first inequality follows from the choice ot Xf 
and P according to Assumption |T] and the fact that 
Bk{x*x) < Bk{,x*x) < x*xMx*x Vx^ G Xf. The sec¬ 
ond inequality follows from the monotonicity of the re¬ 
laxed logarithmic barrier function Bfi-) and the assump¬ 
tion that ip{Akx*^) < <p(x)^) 'ix*^ G Xf. Based on ar¬ 
guments as in the proof of Theorem (TJ this result can be 
used to prove the decrease property (l 20 l) . 
iii) The fact that the value function decreases shows that 
; d) < Jffixo] d) < and hence Xq G Xn{S) for 
any xq G Xn{S). By repeating this argument, the result¬ 
ing closed-loop system state satisfies x{k) G T/v(d) Vfc > 0 
for any x(0) G A/v(d), which shows that all future states 
and inputs will be strictly feasible. Moreover, due to the 
design of the relaxed barrier functions, J^(x; d) is a well- 
defined, positive definite, and radially unbounded func¬ 
tion. Hence, in combination with (uni) it can be used as 
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a Lyapunov function, proving asymptotic stability of the 
origin with a guaranteed region of attraction of at least 
XNiS). □ 

The following results state some useful properties of the 
region of attraction X]s[{S) and show that the feasible set 
Xjy of the corresponding non-relaxed formulation can be 
recovered by making the relaxation parameter arbitrarily 
small. 

Lemma 2 (cf. [l^b Let the assumptions in Theorem\^ 
hold and let the set X]sf{S) be defined according to Defini- 
tion\^ Then, Xn{S) is a nonempty compact and convex 
set. Furthermore, Tjv((5) C Xfj and X]s[{S) —>■ X^ as 
6^0. 

Proof. For any <5 G satisfying Assumptions [2] and [Sj 
is a positive definite function that is convex and 
radially unbounded. As P{S) is strictly positive (see Defi¬ 
nition [S]) , and e G , the first part follows immediately. 
As shown in the proof of Theorem ^ any initial con¬ 
dition xq G Xpf{6) results in strictly feasible input and 
state sequences, which implies that Xn^S) C Xf^. We 
now show that Xjq{5) will contain any compact subset 
of Xfi if we make the relaxation parameter arbitrarily 
small. Assume ccq G Xf^ and let u*(xo), x*(xo), Jff(xo) 
denote the solution of the corresponding nonrelaxed prob¬ 
lem formulation Q. Then, there always exists (5o(a;o) := 
min{-Cfxkixo) -I- -Clukixo) + dl^, 1 - ip{x%{xo), 
* = j k = 0,...,N-1} G 

with the property that the solutions of relaxed and non¬ 
relaxed formulation will be equivalent for all 6 < (5q(xq), 
Le., u*(xo) = u*(xo), x*(xo) = x*(xo), J^(xo;(5) = 
Jfr(xQ). Note that (5o(2;o) is a continuous function of xq 
since both u* (xq) and x*{xo) are continuous due to the 
smooth problem formulation (TT^ . Let us further define 
5q{xo) := max{d G | <5 < do(a:^o), dj(r(xo) < £^(<5)} 
with the property that xq G Xn{5) Vd < dQ(xo). As 
= Jf^{x^-,5^{x^)) = Jn{xo) for all <5 < (5o(a;o) 
and /3((5) is, on the other hand, continuous and strictly in¬ 
creasing for decreasing S, we can state that dQ(a;o) G R^ 
exists for any Xq G Xfj. Moreover, since both (5o(a;o) 
and Jff{xQ) are continuous, S'q^xq) also is a continuous 
function of xq. Consider now an arbitrary compact set 
Aq ^ Xf^ and define do G as do = min 2 ;ogAro 
Due to the continuity of Sq(xo) and the compactness of 
Xq, this value always exists and it holds that Xq C Xn{S) 
for all d < do- This implies that XNiS) —>■ Xn as d —>■ 0, 
which completes the proof. □ 

Corollary 1. Consider the relaxed barrier function based 
MFC problem formulation il‘A) with F(-) given by D61) . 
For any compact set Xq C Xff there exists a do G R+i- such 
that for any d < do and any a;(0) G Xq the feedback uik) = 
UQ{xik)) asymptotically stabilizes the origin of system (jf) 
under strict satisfaction of all input and state constraints. 

It has to be noted that the above conditions for closed- 
loop stability and constraint satisfaction are of course only 


sufficient and may be rather conservative. In particu¬ 
lar, for a given d G R-|+ the actual region of attraction of 
the closed-loop system may be considerably larger than 
the set A/v(d). Likewise, a very small d may be needed 
to achieve Xq C XNiS) when Xq approaches X^. How¬ 
ever, despite possible practical limitations, the presented 
results provide interesting insights and a theoretical jus¬ 
tification for the use of relaxed barrier functions in the 
context of MFC. 

B2) Global stabilization with a nonrelaxed terminal set 

Assuming controllability of system o, we will in the fol¬ 
lowing present a second approach that allows to guarantee 
asymptotic stability for any initial condition by relaxing 
the barrier functions of state and input constraints while 
strictly enforcing the terminal set constraint. 

Assumption 4. The pair (A, B) is controllable and the 
prediction horizon satisfies the condition N > n. More¬ 
over, the terminal cost function is given as 

Fix) = x^PxeBfix), ( 22 ) 

where P and Bfi-) are chosen according to Assumptionl^ 

Theorem 3. Let Assumptions Ql and hold true 
and consider problem m with the terminal cost F(-) 
from i22\) . Then, independently of the relaxation parame¬ 
ter S G R+I-, the feedback uik) = Mo(a:(fc)) asymptotically 
stabilizes the origin of system m for any initial condition 
a;(0) G R". 

Proof. Due to the above controllability assumption, there 
exists for any xq G R” an input sequence ■u(a;o) = 
{mo, ..., UAT-i} such that XNiuixQ), xq) G X^. Hence, 
u*ixo) and J)(,(a;o;d) are defined for any xq G R". 
Since the corresponding terminal state satisfies x*f^ixo) G 
Ay, the local controller u = Kx can again be used 
to construct a feasible control sequence u'^ixo) = 
{mj, ..., u*j,f_i,Kx*j,f} for the successor state a:)}' = Axq -\- 
Buq. Since all parameters are chosen according to As- 
sumption[T] the same arguments as in part ii) of the proof 
of Theorem can be used in order to show that 

J^(4;d)-J^(xo;d) < -ioixQ,u*ixQ))yxo G R" . (23) 

Moreover, due to the design of all involved barrier func¬ 
tions, Jffix;S) is a well-defined, positive definite, and ra¬ 
dially unbounded function. Hence, Jffi-;S) can be em¬ 
ployed as a Lyapunov function for proving global asymp¬ 
totic stability of the origin, see part iii) of the proof of 
Theorem O □ 

The above result illustrates how we can achieve stabi¬ 
lization of the origin for any initial condition, which makes 
the corresponding MFC scheme very robust against un¬ 
certainties, disturbances, or even infeasible start configu¬ 
rations. As we need to strictly enforce the terminal set 
constraint, the benefit of global stabilization comes with 


the cost of possible violations of the state and input con¬ 
straints. However, using the arguments of the previous 
section, we can for any Xq C X'^ still recover strict sat¬ 
isfaction of all state and input constraints by making the 
relaxation parameter <5 arbitrarily small. 

C. Closed-Loop Stability: Terminal Set Free Approaches 

In the previous section, we have seen how asymptotic sta¬ 
bility and even strict constraint satisfaction can be guar¬ 
anteed in the context of relaxed logarithmic barrier func¬ 
tion based MFC by making use of suitable terminal set 
formulations. As in conventional MFC schemes, the ter¬ 
minal set is on the one hand used for ensuring the ex¬ 
istence of a feasible local control law and, thus, the re¬ 
cursive feasibility of the corresponding open-loop optimal 
control problem. On the other hand, only restricting the 
terminal state to a compact set around the origin allowed 
us to derive a quadratic upper bound for the barrier func¬ 
tion see Assumption [T] and the proof of Theorem[I] 

In the following, we will show that the use of relaxed log¬ 
arithmic barrier functions in fact allows us to circumvent 
these two problems and to design novel MFC approaches 
which not only eliminate the need for an explicit terminal 
set constraint but also to prove global asymptotic stabil¬ 
ity of the origin. 

Cl) Tail sequence based terminal cost function 

It is a well-known result that closed-loop stability of both 
linear and nonlinear MFC schemes may be ensured by 
choosing the terminal cost as a suitable CLF that is an up- 
per bound for the infinite-horizon cost-to-go, see e.g. [iq . 
In the presence of input and state constraints, deriving 
such a function in global form is generally not possible, 
which directly motivates the use of a local CLF in com¬ 
bination with a corresponding terminal set constraint. 
However, when considering the relaxed problem formu¬ 
lation (nm, any input sequence which steers the state to 
the origin in a finite number of steps can be used to derive 
an upper bound on the infinite-horizon cost-to-go. 

Assumption 5. Let v{x) := {z;o(a::), ■ • ■ ,vt-i{x)} be an 
input sequence which steers the state of system m to the 
origin in a finite number of T > n steps for any x G R." 
and assume that vi{x) = 0 V/ = 0, ...,T — I x = 
0. Furthermore, let z{x) := {zo(x),..., zr-iix)} with 
zo{x) = x, zi+i(x) = Azi{x) -F Bvi{x), I = 0,...,T -I, 
and zt(x) = 0 be the corresponding state sequence. 

Based on Assumption [5] we propose to choose the ter- 
minal cost as 

T-l 

F{x) = J2kzi{x),vi{x)), (24) 

where £ : R" x R™ —>■ R_|_ refers to the already introduced 
modified stage cost based on relaxed logarithmic barrier 
functions for the state and input constraints, cf. (TT^ . Due 


to the relaxation, v(x) is a valid but possibly suboptimal 
input sequence that steers the state from a;(0) = a; to the 
origin in a finite number of steps. Consequently, F{-) is 
an upper bound for the infinite-horizon cost-to-go, which 
allows us to state the following stability result. 

Theorem 4. Let the Assumptions\^and\^hold and con¬ 
sider problem ilS\} with the terminal cost F{-) from |^]). 
Then, independently of the relaxation parameter S G R^; 
the feedback u{k) = ic^{x{k)) asymptotically stabilizes the 
origin of system m for any initial condition x(0) € R”. 

Proof. Due to Assumption [51 there exist for any Xq G R" 
and any input sequence u(xo) with resulting terminal 
state xn = xn(u(xo),Xq) suitable tail sequences v{xn) 
and z{xm). Moreover, as all barrier functions for the 
state and input constraints are relaxed, see Assumption|2l 
both the stage and terminal cost are always well defined. 
Hence, also u*{xq) and Jff{xo;S) are defined for any 
Xq G R”, which shows that problem (fT^ always admits a 
feasible solution. 

We now want to show that the value function Jff{x{k); 6) 
decreases under the applied feedback for all a;(0) G 
R". Let u*{xo) = {tig,..., and x*{xo) = 

(xo, xl,..., x^} denote the optimal open-loop input and 
state sequences for a given initial condition a;o G R" 
and let v*{xo) := v{x^{xo)) and z*{xo) := z(x’^(xo)) 
be the corresponding tail sequences for the resulting pre¬ 
dicted terminal state. Consider now the successor state 
x'^ =Axo + Bul. Clearly, uA^Xq) = {m{, ..., ii5^r_p-Co} 
and {,..., 0} are suboptimal input se¬ 

quences which steer the state from x^ to the origin in 
a finite number of A -|- T — 1 steps. The resulting state 
sequences are given by a;+(a;o) = {x*,... ,xA, zf} and 
z'^{xo) = {zl,... ,zf<_i,Q}. Using the above subop¬ 
timal input and state sequences and the fact that the 
tail sequences within the terminal cost are simply ap¬ 
pended by zero values, it is straightforward to show that 
Jn{u'^{xo),x'^-,5)-J%{xq-, 5) = -£{xQ,ul{xQ))'ixQ G R". 
From the suboptimality of 'uA{xq) it follows immediately 
that 

J^{4;(5)-J)V(xo;<5)<-£(a:o,fiS(^o))VxoeR" . (25) 

By the design of the relaxed barriers and the assumption 
that vi{x) = 0 'iI = 0,... ,T — 1 FF X = Q, the function 
is well-defined, positive definite, and radially un¬ 
bounded. Thus, it can be used as a Lyapunov function 
for the closed-loop system, proving global asymptotic sta¬ 
bility of the origin. □ 

Note that in order to ensure convexity of the termi¬ 
nal cost function F{-), the elements of the parametrized 
tail input sequence v{x) should be affine in the argu¬ 
ment X. In combination with the condition that vi{x) = 
OVZ = I,...,T— I X = 0, this in fact limits v{-) 
to contain a sequence of linear state feedback laws, i.e. 
v{x) = {Kqx, ... ,Kt-ix}. In the following, we briefly 
discuss two different design approaches for v{x) that meet 
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this requirement and thus allow to guarantee stability of 
the closed-loop system as well as convexity of the result¬ 
ing overall cost function. 

The easiest way to design suitable tail sequences «(•) 
and z(-) is by making use of a linear dead-beat con¬ 
troller. To this end, we may choose a terminal control gain 
K € in such a way that the matrix Ak = A + BK 

is nilpotent, i.e., that it satisfies A^j^ = 0 for some r < n. 
Under the assumption of controllability, this can for ex¬ 
ample be achieved by a suitable pole placement procedure 
which ensures that all eigenvalues of the matrix Ak are 
located at the origin of the complex plane. Based on 
these ideas, we may set T = n and choose the tail input 
sequence as 

v{x) = \^Kx,...,K{A + BKf~^x^ (26) 

which results in the corresponding state sequence z{x) = 
{x,..., A^^x}. Due to the design of the matrix K, it ob¬ 
viously holds that Zn = A^x = 0. Thus, v{-) and z{-) sat¬ 
isfy the conditions in Assumption [SJ and Theorem 2] can 
be used to conclude stability of the closed-loop system. 
Note that, depending on the algorithm, the pole place¬ 
ment problem may become numerically ill-conditioned 
when assigning all poles to exactly the same location. 
Hence, for practical implementations it may be meaning¬ 
ful to distribute the poles of Ak in an e-ball around the 
origin. 

While the above approach allows for a rather simple de¬ 
sign of the tail sequences v{-) and z{-), and thus of the 
terminal cost F(-), the implicit requirement that the pre¬ 
dicted terminal state is steered to the origin in at most n 
steps might be restrictive, leading to suboptimal or even 
aggressive behavior of the overall closed-loop system. In 
the following, we present a second approach that elimi¬ 
nates this restriction by allowing for tail sequences with 
T > n elements. In order to enforce in addition a certain 
optimality with respect to the underlying performance 
criterion, we propose to choose the parametrized tail in¬ 
put sequence v{-) as the solution to the finite-horizon 
LQR problem with zero terminal state constraint, i.e., 

T-l 

v{x) = argmin ^{zi,vi) (27a) 

s. t. Zi+i = Azi + Bvi, I .. ,T - I, (27b) 

zt = 0, zq = X (27c) 

for the standard quadratic stage cost £{z,v) = \\z\\q + 
||u||fj. It can be shown that the solution v*{x) to prob¬ 
lem can be expressed as a sequence of static lin¬ 
ear state feedbacks of the form v*{x) = Kix with I = 
0,...,T — 1. Explicit expressions for the optimal con¬ 
trol law and the corresponding state (and costate) tra¬ 
jectories were given in |26l | based on the solution of two 
unconstrained infinite-horizon LQR problems and a suit¬ 
able iteration scheme. Furthermore, in 0 the authors 


present explicit solutions for more general start and end 
point constraints based on a parametrization of all the 
solutions of the extended symplectic system. However, 
as discussed in the Appendix, the optimal input sequence 
can also be computed directly in vector form as V*{x) = 

[uS^(x) • • • 4li(a;)]^ = Kvx with Ky S In 

both cases, the resulting terminal cost function can be 
formulated as 

T-l 

F{x) = x^Px + e B^{zi{x)) -I- Bu{vi{x)) , (28) 

1=0 

where zi{x) and vi{x) denote the elements of the re¬ 
spective state and input tail sequences and the matrix 
P S can be constructed by inserting these sequences 
into the quadratic stage cost ({■,■), see Appendix. As 
for the previous approach, global asymptotic stability of 
the closed-loop system can be concluded from TheoremSl 
However, the possibility to choose the length of the tail 
sequences within the construction of the terminal cost 
typically leads to an overall improved closed-loop perfor¬ 
mance. Especially if no constraints are active or violated 
by the open-loop state and input tail sequences, the ter¬ 
minal cost function F{-) based on (l27ll may give a quite 
good approximation of the real infinite-horizon cost-to-go. 

Remark 2. Note that the diseussed approaches force the 
predicted state to the origin in N + T steps and are hence 
similar to MPC approaches that are based on an explicit 
zero terminal state constraint, see e.g. Ml- However, by 
using the proposed terminal cost function, this behavior 
is enforced implicitly, i.e., without introducing an explicit 
equality constraint within the optimization problem. Fur¬ 
thermore, stability can be guaranteed for a global region of 
attraction and the tail sequence horizon T may be made 
arbitrary large without increasing the number of optimiza¬ 
tion variables. 

C2) Quadratic terminal cost 

In the previous section, we exploited the fact that the re¬ 
laxation of input and state constraint allows to apply any 
sequence of inputs at the end of the prediction horizon. 
In the following, we will show that in the presence of re¬ 
laxed barrier functions we can in addition derive a global 
quadratic upper bound for the combined state and input 
constraint barrier function, which makes it possible to use 
a purely quadratic terminal cost function term without 
the need for a corresponding terminal set constraint. 

Assumption 6. Let Bxi') and Hu(-) be relaxed gradient 
or weight recentered logarithmic barrier functions accord¬ 
ing to Assumption m see hlS\) and respectively. Let 
the relaxing function be quadratic and given by P 2 {-\d) 
from nUi) . 

Lemma 3. Let {A,B) be stabilizable and let K G 
be a corresponding stabilizing linear control gain. Fur¬ 
thermore, let Assumption\^ hold and consider Bk{x) = 
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B^(x) + B^(Kx) for a given 6 € K 4 f • Then, it holds that 

Bk{x) < x'^{M,, + K^M,,K)x VxGM”, (29) 

where the matrices Mx £ and G j^mxm 

defined as Mx := ^Cjdia.g{t + wfij and M^ := 
^Cj^diag (1 + Wu) Cjj, respectively. Here, Wy, G and 
Wu G R«“ are suitable weighting vectors when consider¬ 
ing weight recentered barrier functions, whereas Wy = 0, 
Wu = 0 in the gradient recentering case. 

Proof. We exemplary consider the state constraints and 
show that Bx{x) < Vx G M". Based on Tay¬ 

lor’s Theorem (see [2^, Theorem 2.1]) we know that 
Bxix) = Bx{0) [Vi3x(0)]^x -I- ix^V^i3x(Ax)x for some 

A G (0,1). Due the recentering of the barrier func¬ 
tions, the first two terms vanish and we get Bx{x) = 
^x^V‘^Bx{Xx)x for some A G (0,1). In particular, it holds 
that Bx{x) < x^Mx Vx G M" for any M G satisfying 
y^Bx{x) ^ M Vx G R”. When considering the quadratic 
relaxing function l32{-',5) from (fTOll . the Hessian of Bx{-) 
is given by 


v^Mx) 

= C'Jdiag(iAi(x), 

■ ■ ■ ,Dq^{x))Cx, 

(30a) 



-cix + (£>5 


where Dfix) 

— J {-C^x+d],Y 


(30b) 

K <5^ 

—C^x -1- dx < <5 . 


Note that iCx = 0 in the case of gradient recentered barrier 
functions. Since Dfix) < Vx G R", it follows im¬ 

mediately that V‘^Bx{x) < ^Cj^diag (1-b Wx) Cx Vx G 
R". Combining this upper bound on the Hessian 
with the previous arguments, it follows that Bx(x) < 
x^MxX Vx G R” with Mx = ^Cjdiag (1-b Wx) Cx. 
Similarly, it is straightforward to show that Bu{Kx) < 
x^K^MuKx Vx G R" with Mu as defined above. □ 


Let now the controller matrix K and the matrix P G 
be solutions to the modified Riccati equation 

K = -{R-^B^PB-^eMu)~^ B^PA (31a) 

P = Aj^PAj^ + K\R + eMu)K -b Q -b eMx , (31b) 

where Mx and Mu are dehned according to LemmaO We 
then propose to choose the terminal cost function as 

.F(x)=x^Px, (32) 

where P G S!^ is the solution to (|3Tb i. Note that the 
controller gain K is in principle arbitrary. However, the 
above choice results in a minimal value of the terminal 
cost function and ensures that for £ —>■ 0 or in the absence 
of constraints, K and P will reduce to the solution of 
the unconstrained LQR problem. We can now state the 
following stability result. 

Theorem 5. Let Assumption hold true and con¬ 
sider problem ilS\} with the terminal cost function F{-) 


from iS^) . Then, independently of the relaxation parame¬ 
ter S G R 4 f, the feedback u(k) = UQ{x(k)) asymptotically 
stabilizes the origin of system m for any initial condition 
x(0) G R”. 

Proof. Again, both u*{xo) and Jff{xo;6) are defined for 
any xq G R" due to the relaxed state and input con¬ 
straints. Based on the quadratic terminal cost func¬ 
tion above, it is now straightforward to show that the 
value function Jff{x(k); S) decreases under the applied 
MFC feedback for all x{k) G R". For any xq G R" 
there exist optimal input and state sequences tt*(xo) = 
{mq, • ■ • and x*{xo) = {xq, x^,...,xj^} as well 

as the suboptimal input and state sequences ■u“''(xo) = 
u*j.^_j^,Kx%} and a;+(xo) = {x*,..., x^^, AkX%} 
for the successor state Xg = Axq -b Buq with Jf[{xQ ; 6) — 
Jff{xo;6) < Jn{u~''{xo),x'^, 6) - Jff{xo;S). This gives us 

~\\X^kx%\\% — \\x}^\\p ... ( 33 ) 
+ \\x%\\q-\-\\K x%\\'jf^-\-eBK{x%) <0 V xj,r G R " , 

where we used the global quadratic bound on Bk{ ) from 
Lemma [3] and the choice of the terminal cost matrix P, 
see dSIb). Thus, (5) can again be used as a Lyapunov 
function for proving global asymptotic stability. □ 


In summary, we can conclude that the concept of re¬ 
laxed logarithmic barrier functions allows us to design 
globally stabilizing MFC schemes without making use of 
an explicit terminal set or state constraint. Instead, the 
presented approaches are based on a suitable design of 
the respective terminal cost function term that heavily 
exploits the properties and advantages of the underlying 
barrier function relaxation. Thus, the presented results 
may be seen as a barrier function based counterpart to 
existing terminal set free MFC approaches relying on suit¬ 
able terminal cost functions ( jlsL l2Ij ), a sufficiently large 
prediction horizon (i0,|3l), or particular controllability 
assumptions (ji, 15, 3). 

However, with the exception of Theorem no guaran¬ 
tees on the satisfaction of state and input constraints have 
been discussed so far. In fact, the presented global stabil¬ 
ity results can only be achieved since the relaxed barrier 
functions allow for in principle arbitrarily large constraint 
violations. In the next section, we will show that the max¬ 
imal violation of state and input constraints in closed-loop 
operation is always bounded, and we discuss how we can 
compute and control it a priori for a given set of initial 
conditions by adjusting the relaxation parameter 5 G R^- 


D. Maximal Constraint Violation Guarantees 

As outlined above, one of the main advantages of model 
predictive control is given by its ability to deal with in¬ 
put and state constraints. On the other hand, possible 
violations of the corresponding constraints seem to be in¬ 
herent to the concept of relaxed barrier functions. In 
the following, we will show how guarantees on the strict 
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satisfaction or the maximal violation of input and state 
constraints can be given for the globally stabilizing re¬ 
laxed logarithmic barrier function based MPC schemes 
from the previous section. We start with the following 
Lemma which allows to upper bound the values of the 
relaxed barrier functions in closed-loop operation by an 
expression that depends on the initial condition. 

Lemma 4. Let problem be formulated based on 

one of the globally stabilizing MPC schemes discussed 
above (Theorems 0 0). Furthermore, let Xd = 

{a:(0), a;(l),... } and u^i = {m( 0), n(l), •. • } with u{k) = 
u(,{x{k)) for k > 0 denote the resulting state and input 
trajectories of the closed-loop system. Then, for any ini¬ 
tial condition Xg = x(0) € M" and any k >0 it holds that 

-Bx(a:(fc)) < i (^JNixg-, S) - , (34a) 

Bu{,u{k)) < i (^Jff{xo-,S) - xJP*^ xo) , (34b) 

where P*^ G S!)^ is the solution to the discrete-time alge¬ 
braic Riccati eguation related to the infinite-horizon LQR 
problem. 

Proof. Based on the respective Theorems we know that 
Jff{x{k -\- 1);(5) — Jff{x(k); S) < —i{x{k),u{k)) for any 
fc > 0 and any x(0) G K.". In particular, this ensures 
that the closed-loop system is asymptotically stable and 
that limfc_^oo d^(x(k); 6) = 0. Summing up over all future 
sampling instants and using a telescoping sum on the left 
hand side, we get that J)(f(a;(0); ^) > i{x{k), u{k)). 

Furthermore, Y^'^=g K^ik), u(k)) = ^{x{k), u(k)) + 

+ Bu{u{k)) and ^ ^ 

a;(0)^P*(. a;(0) due to the optimality of the unconstrained 
infinite-horizon LQR solution. In combination this yields 

OO 

J*Nixi0):5) > x{0 fP:,x{0)+sY,M^{k))+Mu{k)) 

(35) 

and finally, since all terms in the sum on the right hand 
side are positive definite, eB^{x{k)) < j)^(a;(0); i5) — 
j;(0)^P*(.a;(0) as well as eBu{u{k)) < Jff{x{0);6) — 
a;(0)^P*^a;(0) Vfc > 0. □ 

For ease of notation, let now the scalar Q!(a;o; S) G M+ be 
defined as 

a{xo;S) := Jff{xo;6) - Xq P*^Xo . (36) 

As the relaxed barrier function are positive definite and 
radially unbounded. Lemma |4] implies that also the vio¬ 
lations of the corresponding constraints are bounded. In 
particular, for any e G K 4 f, S G K 4 f and any initial con¬ 
dition xq = x(0) G K", upper bounds for the maximal 
violations of state and input constraints are given by 

i*(a;o;(5) =max|(7*x - I Bx(x) < | ( 37 a) 

z!({xo]S) =max I i3u(M) < j. ( 37 b) 


for f = 1,..., gx and j = 1,... ,qu, where d(a;o; d) is de¬ 
fined according to (1551) . Based on this, we can formulate 
the following result on the maximal constraint violations 
that can occur in closed-loop operation. 

Theorem 6. Let problem il2\} be formulated based on one 
of the globally stabilizing MPC schemes discussed above 
(Theorems 0 0) and let Xd = {a;(0), a;(l),... } and 

Ud = {■w(O), m(1), ... } with u{k) = UQ{x{k)) for k > 0 
denote the resulting closed-loop state and input trajecto¬ 
ries. Then, for any xg = a;(0) G K” and any k > 0 it 
holds that 


I^x^(^) ^ dx-(--2^x(^0;d), C,]^u{lf) ^ du-(--2u(^0;d), (38) 

where the elements of the maximal constraint violation 
vectors Zx{xg;6) G K.®* and Zu(a^o;d) G K*" are given 
by 

Proof. Follows directly from Lemma 0] and the definition 
of d(xo;d) and z^(xo;d), z({xg-,6) in (1551) and (1571) . □ 


Note that the optimization problems in (I37|) are con¬ 
vex as Bx(‘) and dlu(’) are convex functions. Moreover, 
when considering gradient recentered logarithmic barrier 
functions, we can exploit the fact that Bx(-) and Bu(-) 
consist only of positive definite terms. In this case, upper 
bounds for the maximal constraint violations are given by 
z^{xg-,S) < max{—Zx,0}, where z^. is given as 


= mm 


z\l3{z-,5) +\n{d]^ 


- 1 = 


a{xg-,5) 


(39) 

for i = 1,..., gx and similar for the input con¬ 

straints, cf. (1T51) . If the relaxing function /3(-; 6) is chosen 
as (3e{-',5) or /3fe(-;d) with k > 2, a. nonlinear equation 
solver can be used to find suitable solutions z®. However, 


in the case of /3(-; d) = /3fe(-; d) with k = 2, the equality 
constraint in (1551) reduces to a quadratic equation in z 
and a closed form expression of the maximal constraint 
violation can be given as 


4 = - \Jlli - 7*, 2 ) , f = 1,..., gx , (40) 

where 7 i,i :=2-jr and 7 i ,2 := 1-1-2 In - |d(a;o;d), 
see also Lemma 21 in [l^. Note that the values for the 
maximal constraint violations obtained above may also be 
negative. In this case, the corresponding constraints will 
not be violated but rather satisfied with the respective 
safety margin. An important consequence of Theorem [5] 
is that, as stated in the following Lemma, there always 
exists a set of initial conditions for which the input and 
state constraints will not be violated at all. 
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Lemma 5. Let problem U2\) be formulated based on one 
of the globally stabilizing MPC schemes discussed above 
(Theorems and let Xd = {a;(0), a;(l),... } and 






Uci = {m( 0), u(l),... } with u{k) = UQ{x{k)) for k > 0 de¬ 
note the resulting closed-loop state and input trajectories. 
Furthermore, for S G let the set be defined as 

:= {x G X \ a{x, S) < £^'(5)} , (41) 

with P'{S) := min{^x(<5), ,Su(^)}; where Pxid) and ^u(^) 
are defined according to Definitions^ Then, for any ini¬ 
tial condition a;(0) G Xfi{5) and any k > 0 it holds that 
Cxx{k) < dx as well as Cuu{k) < du- 

Proof. For any xg = x(0) € and all fc > 0 it 

holds due to Lemma 0] and the definition of f'{S) that 
eBy,[x{k)) < d(xo;(5) < and eB,i{u{k)) < a{x,S) < 

e/3u(d). In combination with the definition of /3x(^) and 
^u(^), it follows directly that x{k) € X, u{k) G U Vk > 0, 
cf. Lemma [TJ □ 

One may now ask how large we can make the set of 

initial conditions for which strict satisfaction of all input 
and state constraints is guaranteed. In the following, we 
give an answer to this question for each of the different 
MFC approaches discussed above. In particular, we show 
that it is always possible to recover the feasible set of a 
suitable corresponding nonrelaxed MFC formulation. 

Theorem 7. Let problem S12\) be formulated based on one 
of the globally stabilizing MFC schemes discussed above 
( Theorems [3 0 0^- Moreover, let Xm be defined as 

Xn = {x G X :3u s.t. Xk G X,Uk G U,xn S Xf} (42) 

when considering the approach based on a nonrelaxed ter¬ 
minal set constraint (Theorem\B^, as 


When considering the approach based on Theorem [5] the 
problem occurs that the matrix P from (13111 . and hence 
the quadratic terminal cost function F{x) = x^Px, will 
grow without bound for d —>■ 0. Also, in this case it 
does not necessarily hold that Xfi^S) C Xn as the re¬ 
sulting optimal state and input sequences for a given 
Xg G Xfj{5) will in general not satisfy the additional 
constraint Xn = 0. However, for any xg G Xf^ there 
exists input and state sequences ■u(xg) and x(xg) which 
strictly satisfy the conditions specified in (I44L in par¬ 
ticular xn(u(xo), Xo) = 0. For these sequences, we can 
always find a dg(xg) G R-|+ such that Jn(u(xo), xg; S) = 
J]v(u(xo), Xg) for all d < do(xg), where J]s[iu{xo),xo) de¬ 
notes the value function of a nonrelaxed problem formu¬ 
lation with the additional constraint xjv = 0, cf. the 
proof of Lemma [2l Furthermore, due to the suboptimal¬ 
ity of the input sequence ■u(xg), it holds that J)(f(xo; d) < 
JN{u{xg),xg-,d). Hence, Jn{xo;S) < JN{uixo),xo) for 
all d < dg(xo), which shows that for any xg G Xf^ the 
value function J)(f(xg;d), and hence also the expression 
d(xg;d), will stay bounded as d —>■ 0. As, on the other 
hand, /3(6) increases without bound and both do(xo) and 
Jn{xo;S) are continuous functions, there exists for any 
compact set Xg C Xfj a dg G < min 3 ,gg;t’o ^o(a;g) 
such that Xg C Xf^{6) Vd < dg, see also the proof of 
Lemma 01 □ 

Note that for a general stabilizable system, the set Xn 
in (|4^ may be restricted to a lower-dimensional subspace 
of R” or it may even be empty. However, in case of a 
controllable system it is straightforward to show that Xn 
is a nonempty polytope. 


Xn = {x G X -. 3u s.t. Xk G X ,Uk GU, 

zi{xn) G X,vi{xn) G U} 

when considering the approach based on auxiliary tail se¬ 
quences zif) and vi{-) with I = 1,... ,T — 1 (Theorem^, 
and as 


Xn = {x G X :3u s.t. Xk G X,Uk G U,xn = 0} (44) 

when considering the approach based on a purely quadratic 
terminal cost function (Theorem \^, where in all three 
cases Xg = X and Xk = Xk{u,x) for k = 1,..., N . Then, 
for any compact set of initial conditions Xg C Xff there 
exists a Sg G such that Xg C X(f{S) for all S < Sg. 

Proof. We first consider the approaches based on Theo¬ 
rem |3] and Theorem m The proof is closely related to that 
of Lemma[2]and only a sketch is given here. In particular, 
it can be shown that Xff{d) is a nonempty and compact 
set with Xfj{6) C X^ and Xf^{6) —>■ Xn as d —>■ 0. The 
result is again based an the fact that the existence of a 
strictly feasible solution implies that for any xg G Xfj 
there exists a dg(xg) such that d(xg; d) will stay constant 
for all d < dg, whereas /3(d) can be made arbitrarily large 
as d —y 0. 


In general, the parameter dg that is sufficient for ensuring 
strict constraint satisfaction based on the above results 
may be very small. For practical applications, it might 
therefore be reasonable to enforce the satisfaction of input 
and state constraints only with a predefined tolerance. 
Based on the above arguments, we can easily state the 
following result. 

Corollary 2. Let problem m be formulated based on 
one of the globally stabilizing MPC schemes discussed 
above (Theorems\^\^\^ and let the respective sets Xn be 
given according to Theorem^ Then, for any given con¬ 
straint violation tolerance ztoi G M+ and any given set of 
initial conditions Xg C Xf^ there exists a Sg G M-|+ such 
that for all d < dg and any initial condition x(0) G Xg 
the maximal possible constraint violation is less than ztoi, 
i.e., that for all k > 0 


C^x{k) < dx + ztoi^, C,j,u{k) < du + Ztoi^ ■ (45) 

Based on the above results, we can formulate the follow¬ 
ing iterative algorithm that allows to determine a priori a 
sufficiently small relaxation parameter for ensuring satis¬ 
faction of a given maximal constraint violation tolerance 
for a given set of initial conditions. 
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Algorithm 1 (Compute <5o for Xq and ztoi & R +) 

Input: problem formulation, set Xq, tolerance ztoi 
Output: parameter G s.t. the maximal con¬ 

straint violation is bounded by ztoi for any x( 0 ) S Xq 
1 : choose initial e,S £ R+i- and set up problem (IT^ 

2: repeat 

3: L decrease relaxation parameter: <5 ^ yd, 7 S (0,1) 

4: L determine a{XQ,S) = iiiaxxGXo^ix,S) based 

on ([36)1 

5: L compute zl.{Xo,6), £;j(d:b,d) from (1571) with 

a(A’o,d) 

6: until zl.{Xo;S) < Ztoi and zl{Xo;6) < Ztoi holds V i,j 


Remark 3. Up to now, it has not been clarified whether 
the function a{x\ 6) that is maximized in step 4 of Al¬ 
gorithm 1 is convex. However, for all tested parameter 
configurations and examples with convex Xq , the maximal 
value of a{x; S) was in fact always attained at on of the 
vertices of the set Xq . Moreover, a conservative solution 
can always be found by evaluating a convex upper bound 
ofa{x;S), e.g., the value function Jff{x]6) itself, at the 
vertices of the set Xq . 

Summarizing, we can state that, despite the use of the 
relaxed barrier functions, the maximal possible violation 
of input and state constraints in closed-loop operation is 
bounded and depends directly on the choice of the re¬ 
laxation parameter 5. Furthermore, the presented results 
allow to compute an estimate for the maximal constraint 
violation a priori or to determine a suitable relaxation pa¬ 
rameter that guarantees satisfaction of a given constraint 
violation tolerance. Note that certain constraints can be 
prioritized by making use of different 5i in the relaxation. 
This important fact may for example be used in order 
to enforce satisfaction of physically motivated hard input 
constraints. 

E. Closed-Loop Performance and Robustness 

In this section, we briefly discuss some aspects concern¬ 
ing the performance and robustness properties of the pre¬ 
sented relaxed barrier function based MFC approaches. 
Of course, a thorough investigation of these issues is well 
beyond the scope of this paper and can be considered 
as future work. We begin with the following arguments, 
which show that the presented relaxed barrier function 
based MFC schemes will always recover the closed-loop 
performance of a related MFC scheme based on nonre- 
laxed barrier functions if the relaxation parameter 6 is 
small enough as well as that of a conventional linear MFC 
scheme when the barrier weighting e goes to zero. 

Let the relaxed barrier function based MFC problem m 
be formulated based on one of the stabilizing design ap¬ 
proaches discussed above (Theorem[51[31l31 or[S|), and let 
the set A/v be defined according to (US) for the termi¬ 
nal set based approaches and according to (H51) and (H51) 


for the respective terminal set free approaches. Consider 
Jfi{xQ-, 5) for a given xq G X^ and let Jff{xQ) and Jff{xQ) 
denote to the value functions for the corresponding nonre- 
laxed problem formulation and a conventional MFC for¬ 
mulation, both based on the constraints that define the 
set Xn. By definition, there exists for any xq £ X^ an 
feasible state and input sequences that result in strict 
satisfaction of all inequality constraints in the respec¬ 
tive description of A/v, which implies that there exists a 
(5o(a;o) G K+i- such that Jfj{xo;d) = Jffixo) V<5 < (5o(a;o), 
see the proof of Lemma [5] In particular, for any compact 
set of initial conditions Xq C Xf^ there exists a do G R-h- 
such that for any d < do, any a:(0) G Xq, and any A: > 0 it 
holds that Jfj{x{k)-,5) = Jfj{x{k)). This shows that the 
performance of a corresponding nonrelaxed formulation 
can always be recovered within the interior of the respec¬ 
tive feasible set. Furthermore, it is well known that the 
solution of the nonrelaxed barrier function based prob¬ 
lem converges to the solution of the corresponding 
conventional problem when the barrier function weight¬ 
ing parameter e approaches zero ch. 11]. Thus, for 
arbitrary small e, d G K-h-, the presented relaxed bar¬ 
rier function based MFC schemes recover the closed-loop 
behavior and performance of related conventional MFC 
schemes. The question whether and under which circum¬ 
stances the relaxed barrier function based setup may also 
lead to an overall improved performance of the closed-loop 
system, e.g., a decreased cumulated cost, can be consid¬ 
ered as possible future work. To the experience of the au¬ 
thors, the resulting closed-loop behavior is already very 
good for moderate values of the barrier parameters, i.e. 
e, <5 in the order of 10“^, and typically close to the so¬ 
lution of the original MFC problem for e, d in the order 
of 10-4. 

Concerning the robustness of the closed-loop system, it 
has to be noted that, due to the relaxation of the under¬ 
lying state and input constraints, the resulting overall cost 
function is defined for any x £ M". Hence, the presented 
MFC schemes are robust against arbitrary effects caused 
by disturbances, uncertainties, or measurement errors in 
the sense that the corresponding open-loop optimal con¬ 
trol problem always admits a feasible solution. However, 
can we say something about the qualitative robust stabil¬ 
ity properties of the closed-loop system? For example, we 
might consider the disturbance affected dynamics 

x(k + 1) = Ax(k) -\- BuQ(x{k)) -\- w{k ), (46) 

where w{k) G K" denotes an unknown but bounded addi¬ 
tive disturbance at time instant k > 0. When considering 
the approach based on Theoremj^ any stability or conver¬ 
gence guarantees of the closed-loop system will in general 
be lost as the disturbance might cause the system state to 
leave the set Xm{5). On the other hand, the approaches 
based on Theorems [31 SI and [3] allow to guarantee global 
asymptotic stability of the undisturbed closed-loOT sys¬ 
tem. Interestingly, based on the main result of [ij, this 
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directly implies that system (H51) is integral input-to-state 
stable (iISS) with respect to the disturbance w. As an im¬ 
portant consequence, the state remains bounded and con¬ 
verges asymptotically to the origin for any disturbance se¬ 
quence with bounded energy, see [ij for more details. Fur¬ 
ther investigations, e.g., proving stronger input-to-state 
stability results for system (|i51) . are the subject of possi¬ 
ble future work. 

V. Example and Numerical Aspects 

In this section, we outline a selection of steps which can be 
used to get from a given control problem to a stabilizing 
relaxed logarithmic barrier function based MFC scheme, 
and illustrate both the design and the behavior of the 
closed-loop system by means of a numerical example. In 
addition, we also briefly comment on some interesting nu¬ 
merical aspects of the resulting open-loop optimal control 
problem. 

A. Overall MFC Design Procedure 

Step 1: Basic problem setup. Choose suitable values 
for the problem parameters that are not related to the 
barrier functions, i.e., the weighting matrices Q G S" and 
R S as well as the prediction horizon N S N+. 

Step 2: Relaxed barrier functions. Decide on proce¬ 
dures for relaxing and recentering the logarithmic barrier 
functions and choose suitable (initial) values for the bar¬ 
rier function parameters eGM+i-, (jeK+i-. In general, the 
quadratic relaxation P{-',S) = /32{-',S) seems to work well 
in practice. 

Step 3: Terminal cost function. Choose a suitable 
approach that allows to design the terminal cost F{-) in 
such a way that asymptotic stability of the closed-loop 
system is guaranteed. The different approaches presented 
in this paper are summarized in Table [T] together with ad¬ 
ditional information regarding the necessity of a terminal 
set Xf, the underlying assumptions, and the respective 
region of attraction (ROA). 

Step 4: Parameter tuning. Formulate the open-loop 
optimal control problem (TT^ based on Steps 1-3. The 
relevant problem parameters may be adjusted in order to 
achieve a desired closed-loop performance, e.g. based on 
closed-loop simulations. In addition to the usual MFC 
parameters, the parameters e and S as well as the chosen 
terminal cost function will have a major impact on the 
resulting behavior. To the experience of the authors, and 
in accordance with the above results, the barrier func¬ 
tion weighting parameter e may be used to influence the 
closed-loop performance while the relaxation parameter 6 
primarily allows to control the occurrence and the max¬ 
imal amount of state and input constraint violations. In 
particular. Algorithm 1 from Section [D] may be used to 
adapt the relaxation parameter in such a way that sat¬ 
isfaction of a given constraint violation tolerance can be 
guaranteed. 


B. Numerical Example 

In the following, we briefly illustrate the outlined design 
procedure and the closed-loop behavior of the proposed 
MFC schemes by means of an academic numerical exam¬ 
ple. We consider a discrete-time double integrator system 
of the form 


^{k -I-1) = 


'1 Ts 


'rp2' 

0 1 

x{k) + 

^ S 

T 

s 


i(fc), 


T, = 0.1. (47) 


with the input and state constraint sets U = {u G K. : 
—2 < M < 1} and A = {a: G : —2 < xi < 3, \x 2 \ < 
0 . 8 }. 

Following Step 1 of our design procedure, we choose the 
basic problem parameters to N = 10, Q = diag(l, 0.1), 
and R — 1. In a second step, we decide on using weight 
recentered logarithmic barrier functions with a quadratic 
relaxation based on / 3 (-; 6) = / 32 (-; S), see Definition[T]and 
Eq. (1101) . We fix the barrier function weighting parameter 
to £ = 10“^ and set up the different MFC schemes from 
Table [T] for varying values of the relaxation parameter S. 
For the terminal set based approaches from Section [B] we 
employ a contractive polytopic terminal set as discussed 
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Inspired by our results above, we are in particular inter¬ 
ested in the region of attraction of the locally stabilizing 
approach from Section El as well as in the regions with 
guaranteed strict or approximate constraint satisfaction 
when considering the globally stabilizing approaches from 
the Sections iBb and[C]l /C2. Furthermore, we would like 
to compare the closed-loop behavior of the corresponding 
different MFC schemes and illustrate how the maximal 
constraint violation can be controlled by adjusting the 
relaxation parameter. 

Fig. [2] depicts and compares some of the ^-dependent sets 
that have been discussed in the results above, i.e., regions 
of initial conditions for which we can guarantee properties 
like asymptotic stability or satisfaction of input and state 
constraints. Note that the superscript in each case in¬ 
dicates the respective MFC approach by referring to the 
corresponding theorem number (we use extra indices a 
and b to distinguish between the two tail sequence based 
approaches). While the sets X'^{S) and X'^{5) are given 
as sublevel sets of the value function J^{x\5), the sets 
X^ are sublevel sets of the resulting maximal constraint 
violation based on Theorem[6l both for a given 5 = 10“®. 
In this case, the functions were evaluated over a fine grid 
and Matlab’s contour function was used for plotting. 
Whereas the sets for which strict constraint satisfaction 
can be ensured may be very small, approximate constraint 
satisfaction with a meaningful tolerance of ztoi = 10“^ 
is achieved for much larger regions of initial conditions. 
The approach based on a purely quadratic terminal cost 
F{x) = x^Px results in very small sets which is due to 
the conservative quadratic upper bound (see Lemma El), 
that causes P to grow rather fast depending on the rela¬ 
tion e/(5^. 
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Table 1: Summary of the presented MFC schemes based on relaxed logarithmic barrier functions. 


Approach 

Section 

Terminal cost F{x) 

Xf 

{A,B) 

Relaxing function 

Theorem 

ROA 

constraint violation 

Relaxed Xf 


Px + eBf{x) 

yes 

stabilizable 


Thm. |5] 

Xn(S) 

none 

Nonrelaxed Xf 


x^ Px + eB{(x) 

yes 

controllable 


Thm. 0] 


adjustable, none in Xj^^(6) 

Tail sequences 

n 

YaYo i(zi{x),vi{x)) 

no 

controllable 


Thm.EI 

R" 

adjustable, none in X^{S) 

Quadratic bound 


x^Px 

no 

stabilizable 

02{- ;S), quadratic 

Thm. |S] 

R" 

adjustable, none in X^{S) 


In Fig. 131 the behavior of the resulting closed-loop sys¬ 
tems is illustrated for different initial conditions and a 
varying relaxation parameter 6. It can be seen how con¬ 
vergence to the origin is achieved even for infeasible initial 
conditions and how strict or approximate constraint sat¬ 
isfaction may be enforced by making 6 sufficiently small. 
Note that, in order to get comparable closed-loop perfor¬ 
mance and execution times, we chose a doubled horizon 
of = 20 for the MFC schemes based on a dead-beat 
tail sequence and a purely quadratic terminal cost. 


C. Numerical Aspects 

As outlined above, one of the main advantages of relaxed 
barrier function based MFC formulations is given by the 
fact that the stabilizing control input can be character¬ 
ized as the minimum of a globally defined, continuously 
differentiable, and strongly convex cost function that is 
parametrized by the current system state. In fact, after 
elimination of the linear system dynamics, the open-loop 
optimal control problem can be formulated as the uncon¬ 
strained minimization of a cost function of the form 

.Jn{U, x) = \u^HU + x^FU + x^Yx + eB^u{U, x ), (48) 


where U := ••• S R"^, nu = Nm, and 

Bxu : X R” — ^ R_|. is a positive definite, con¬ 

vex, and continuously differentiable relaxed logarithmic 
barrier function for polytopic constraints of the form 
GU < w + Ex. The matrices i? S F e 
F e S!;:, C? e w e R«, and E e R?^" can be 

constructed from lEl) and the corresponding constraints 
by means of simple matrix operations. Note that the 
matrices in the above QF formulation may be rather ill- 
conditioned when considering unstable systems and long 
prediction horizons. This issue can, however, be resolved 
by a suitable prestabilization, see [3l| . The then resulting 
mixed input and state constraints are not considered here 
in detail but can be handled by all discussed approaches 
after some minor modifications. While the above con¬ 
densed QF formulation allows for a compact represen¬ 
tation with a minimal number of optimization variables, 
the original formulation (jT^ might allow to exploit the in¬ 
herent problem structure, which may be beneficial from a 
numerical point of view. Tailored structure exploiting op¬ 
timization algorithms have been presented in the context 
of both conventional linear MFC, e.g.ls^ and 33|, and 
nonrelaxed barrier function based linear MFC, see 341. 



Figure 2: The discussed sets for the different MFC schemes 
with e = 10 “^ and S = 10 “®. Here, Xf and Xn denote 
the polytopic terminal set and the resulting feasible set of the 
terminal set based approaches. Furthermore, Ty((5) denotes 
the region of attraction of the locally stabilizing approach with 
strict constraint satisfaction, while X'^{5) and denote 

the sets with guaranteed zero constraint violation for the glob¬ 
ally stabilizing approaches based on Theorems\^ and^ where 
the LQR based tail sequence is used in the latter case. The 
corresponding sets X^°'{5) for a dead-beat controller based tail- 
sequence and X§{5) for a purely quadratic terminal cost where 
in this case too small to be plotted. In addition, the much larger 
sets Xq{S), A’o“(5), Xq^{S), and Xq{S) denote the regions in 
which the respective MFC schemes result in a maximal con¬ 
straint violation of Ztoi = 10~®. 


Whether such techniques could also be applied for the re¬ 
laxed barrier function based MFC approaches discussed 
in this paper, is an interesting open question which might 
be considered as possible future work. 

Another interesting property of the condensed formula¬ 
tion (|48l) is that it can be tackled directly by Newton- 
based continuous-time optimization algorithms as they 
are, for example, discussed in [l^ and 121. In partic¬ 
ular, the differentiability and convexity properties of the 
cost function allow to design a continuous-time dynamical 
system of the form 


U{t)= f{Uit),x{t)), U{to) = Uo 


(49) 


whose solution asymptotically tracks the optimal input 
vector U*{x{t)), where x{t) is the continuously measured 
system state, see the aforementioned references for more 
details. Such continuous-time linear MFC algorithms 
completely eliminate the need for an iterative on-line op- 
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Figure 3: Left: closed-loop behavior for * 0,1 = [2.15, —0.7]^ £ Xn and 3 £ {0.1, 2.5 x 10“^, 10~^, 10“^, 10~^} for the approach 
based on a nonrelaxed terminal set [*]: if 3 is small enough, the state trajectories stay feasible. Also shown is the behavior for two 
infeasible initial conditions xo ,2 and 2:0,3 with 3 — 10“® together with the corresponding predicted terminal states x%{x{k)) [of. 
For the infeasible initial condition xo ,2 = [—1.75, —1]^, also the closed-loop trajectories of the globally stabilizing terminal set 
free MFC schemes from SectionY^ are depicted (dead-beat based tail sequence [o], LQR based tail sequence /+/, purely quadratic 
terminal cost [a]). 

Middle: control input for the different globally stabilizing MFC schemes for the initial condition xo ,2 = [— 1 . 75 , — 1 ]^ using the 
same coloring scheme. 

Right: The relaxation parameter Jo obtained by Algorithm 1 for ensuring strict constraint satisfaction [of, i.e. Ztoi = 0, respec¬ 
tively approximate constraint satisfaction with a tolerance of ztoi = 10“® [*] for the approach based on Theorem\^with initial 
conditions in Xq = kXn where k £ (0, 1 ). Also depicted is the So required for ensuring asymptotic stability and strict constraint 
satisfaction directly based on the approach with a relaxed terminal set [o], cf. Theorem\^ 


timization and can be implemented for in principle ar¬ 
bitrary fast system dynamics. We expect that it might 
be possible to exploit some of the advantages of relaxed 
logarithmic barrier function based formulations, e.g., the 
global definition of the cost function and a bounded cur¬ 
vature, explicitly within the underlying numerical inte¬ 
gration. 

VI. Conclusion 

In summary, our investigation showed that the concept of 
relaxed logarithmic barrier function based model predic¬ 
tive control is interesting both from a system theoretical 
and a practical point of view. In particular, while we 
are still able to recover many of the theoretical proper¬ 
ties of conventional or nonrelaxed barrier function based 
MFC schemes, the use of relaxed logarithmic barrier func¬ 
tions allows to characterize the stabilizing control input as 
the minimizer of a globally defined, continuously differen¬ 
tiable and strongly convex function that is parametrized 
by the current system state. As a main result, we pre¬ 
sented different constructive MFC design approaches that 
guarantee global asymptotic stability and allow to influ¬ 
ence the performance and constraint satisfaction proper¬ 
ties of the closed-loop system directly by adjusting the 
relaxation parameter. The resulting MFC schemes are 
not necessarily based on the construction of a suitable 
terminal set and, due to the underlying relaxation, inher¬ 
ently robust against disturbances, uncertainties, or sensor 
faults. 


servative constraint violation bounds, or the design and 
comparison of tailored iterative or continuous-time opti¬ 
mization algorithms that explicitly exploit the properties 
of the relaxed MFC problem formulation. 

Appendix 

In the following, we show how the finite-horizon LQR 
problem with zero terminal state constraint can be solved 
without the use of algebraic Riccati equations. Consider 
problem (12711 for a given initial state x and horizon T > n. 
By eliminating the predicted states zi for I = 1..., T 
using (I27b i. we get the following equivalent formulation 
in vectorized form: 

Jy[x) = min -V^HV -\- -x^FV -\- -V~^F^x -\- 
' V 2 2 2 2 

s. t. A^x + [^i S' 2 ] fo = 0 , (50) 

where the respective matrices are given as 
Si = [A^-^B ■ ■ ■ A"R] , 5*2 = ■■■ B] 

H = 2[r + , F = 2L2^Q$, Y = 2{q + 


■ A' 


- B 

• O' 

Q = It 'S> Q 

A^ 

, $ = 

It- 

4 

1 • • • 

to 

• B 

R = It ® R 


Due to the controllability of the system, there always ex¬ 
ists a feasible solution V{x) to problem (l50t . We parti¬ 
tion the vector V £ as with Vi £ 

]^(T-n)m g K"™. The terminal state constraint 

can be eliminated by choosing V 2 = — 5^ (A'^x -|- S'lVi), 


Interesting open problems may include a thorough anal¬ 
ysis of the performance and robustness properties of the 
resulting closed-loop system, the derivation of less con¬ 
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where S 2 denotes the Moore-Penrose pseudoinverse of the 
matrix 5*2. This leads to 


V = 


^{T—n)m 

-s+s. 


Pi 


0 

-S+A^ 


X = TvVi+T^x . (51) 


Inserting (I5T|) into problem (I50p results in the following 
unconstrained optimization problem in the variable Vi 

Jv{x) = min ^-V^HVi + x^FVi + \-x^Yx , (52) 

V\ z z 


where the transformed problem matrices are given by 
H = Vj,HTv, F = VlHVv + and f = P + 

T^FlTx + + FT^. Note that since V{x) = T^x 

is a feasible solution and Im(ry) = Null([5'i <S'2]), the 
solution is invariant under this change of coordinates, 
and the problems dsni) and ([52l) are equivalent, cf. [^, 
p. 132]. Furthermore, it always holds that H G sfj 
as the matrix Fy has full column rank. The optimal so¬ 
lution to the reduced problem can be computed easily as 
Vi{x) = —H~^F^x which results in 

V*ix) = (r^-rvH~^F^^x = Kvx . (53) 

Furthermore, the optimal cost is given by Jy{x) = x^Pyx 
with Py = \(^ G- F^H-^F^ G S!^. 
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